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Abstract. 

-^ The correlation coefficients of the linear regression of six solar indices ver- 

^—1 sus FiQj were analyzed in solar cycles 21, 22 and 23. We also analyzed the 

r^ interconnection between these indices and -Fio,? with help of the approxima- 

^ tion by the polynomials of second order. The indices we've studied in this 

^ paper are: Wolf numbers - W, 530,3 nm coronal line flux - -F530, the total so- 

"*^ lar irradiance - TSI, Mg II UV-index 280 nm core-to-wing ratio. Flare Index 

^ - FI and Counts of flares. In the most cases the regressions of these solar 

indices versus -Fio,? are close to linear except the moments of time near to the 

minimums and maximums of 11-year activity. For the linear regressions we 

found that the values of correlation coefficients Kcorr{t) for the indices versus 

Fioj and W show the cyclic variations with periods approximately equal to 






P_i the to half length of 11-year cycle - 5,5 years approximately. 
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L^ 1 Introduction 



►>. Magnetic activity of the Sun is called the complex of electromagnetic and 

l/~l hydrodynamic processes in the solar atmosphere. The analysis of active re- 

^ gions (plages and spots in the photosphere, flocculae in the chromosphere 

Tj" and prominences in the corona of the Sun) is required to study the mag- 

■^ netic fleld of the Sun and the physics of magnetic activity. This task is of 

^ fundamental importance for astrophysics of the Sun and the stars. Its ap- 

,-H plied meaning is connected with the influence of solar active processes on the 

j>! Earth's magnetic fleld. 

IT^ We have studied monthly averaged values of six global solar activity in- 

j[^ dices in magnetic activity cycles 21, 22 and 23. Most of these observed data 

we used in our paper were published in Solar-Geophysical Data Reports 
(2009) and National Geophysical Data Center Solar and Terrestrial Physics 
(2012). 

All the indices studied in this paper are very important not only for anal- 
ysis of solar radiation formed on the different altitudes of solar atmosphere, 
but for solar-terrestrial relationships as the key factors of the solar radiation 



C^ 



influence (EUV/UV solar radiation is the most important) on the different 
layers of terrestrial atmosphere also. 

It's known that the various solar activity parameters correlate quite well 
with the most popular solar index sunspot numbers W and with each others 
over long time scales. Floyd et al. (2005) showed that the mutual relation 
between sunspot numbers and three solar UV/EUV indices, the -Fioj flux, 
the Mgll core-to-wing ratio, which is the important chromospheric UV index 
(Viereck et al. (2001), Viereck et al. (2004)), remained stable for 25 years 
until 2000. At the end of 2001 these mutual relations dramatically changed 
due to a large enhancement which took place after actual sunspot maximum 
of the cycle 23 and the subsequent relative quietness intermediate called the 
Gnevyshev gap. 

In our issue for all the indices we used the monthly averages values. Such 
averages allowed us to take into consideration the fact that the major mod- 
ulation of solar indexes contains a periodicity of about 27-28 days (corre- 
sponding to the mean solar rotation period). So we reduced the influence of 
the rotational modulation of the data sets. 

Vitinsky et al. (1986) analyzed solar cycles 18 - 20 and pointed out that 
correlation for spot numbers versus radio flux Fiqj does not show the close 
linear connection during all the activity cycle. Also it was emphasized the 
importance of statistical study in our solar activity processes understanding. 
To achieve the best agreement in approximation of spot numbers values by 
FiQj observations Vitinsky et al. (1986) proposed to approximate the de- 
pendence W - -Fio,7 by two linear regressions: the first one - for the low solar 
activity (where -Fioj less than 150) and the second one - for the high activity 
(Fio,7 more than 150). 

In this paper we found out that the linear correlation was violated not 
only for maximums of solar activity cycles but for minimums of the cycles 
too. Our analysis of the interconnection between these indices and -Fio,? with 
help of the approximation by the polynomials of second order confirmed this 
fact. 

We also analyzed a three- year determined correlation coefficients Kcorr{t) 
for six solar activity indices versus Fiqj. 

2 Global activity indices 

Then we have to say a few words about solar indices studied in this paper. 

At Figure 1 we show the activity indices which are normalized to their 
values averaged over the analyzed time interval. We can see that the relative 
variation of the index in a solar cycle is about 2-3 times. However, the 
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Figure 1: The time series of monthly average values of W, -Fioj, Mgll core- 
to-wing ratio, -Fsao^s, Flare Index and Counts of flares. The upper index N 
indicates that solar activity indices are normalized to their values averaged 
over the analyzed time interval. 



magnitude of Mgll - index (as well as TSI , not shown in this figure ) change 
very little - about shares of percent. 

The solar radio microwave flux at wavelengths 10,7 cm Fiqj has the long 
running series of observations started in 1947 in Ottawa, Canada and main- 
tained to this day at Penticton site in British Columbia. This radio emission 
comes from high part of the chromosphere and low part of the corona. -Fio,7 
radio flux has two different sources: thermal bremsstrahlung (due to elec- 
trons radiating when changing direction by being deflected by other charged 
participles) and gyro-radiation (due to electrons radiating when changing di- 
rection by gyrating around magnetic fields lines). These mechanisms give rise 
to enhanced radiation when the temperature, density and magnetic fields are 
enhanced. So -Fio,? is a good measure of general solar activity. Fiqj data are 
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Figure 2: The time series of the average annual W from 1700 for 2009, and 
its wavelet image (we use the Morley wavelet transfer function). 



available at |http://radbelts. gsfc.nasa.gov/RB_model_model_int/Psi_database.htnil 



The Mgll 280 nm is important solar activity indicator of radiation, de- 
rived from daily solar observations of the core-to-wing ratio of the Mgll 
doublet at 279,9 nm provides a good measure of the solar UV variability 
and can be used as a reliable proxy to model extreme UV (EUV) variability 
during the solar cycle Scupin et al. (2005). The Mgll observation data were 
obtained from several satellite's (NOAA, ENVISAT) instruments. NOAA 
started in 1978 (during the 21"**, 22"*^ and the first part of the 23''"^ solar 
activity cycles), ENVISAT was launched on 2002 (last part of the 23*'' solar 
activity cycle). Comparison of the NOAA and ENVISAT Mgll index obser- 
vation data shows that both the Mgll indexes agree to within about 0,5%. 
We used both the NOAA and ENVISAT Mgll index observed data from 
Puga and Viereck (2004), Scupin et al. (2005). Viereck et al. (2001) showed 
an extremely good fit between the 30,4 nm emission (the main component 



of EUV-emission) and the NOAA Mgll index. Viereck et al. (2001), Floyd 
et al. (2005) have noted the linear relationship between the Mgll index and 
total solar irradiance. 

Solar irradiance is the total amount of solar energy at a given wavelength 
received at the top of the earth's atmosphere per unit time. When integrated 
over all wavelengths, this quantity is called the total solar irradiance (TSI) 
previously known as the solar constant. Regular monitoring of TSI has been 
carried out since 1978. From 1985 the total solar irradiance was observed 
by Earth Radiation Budget Satellite (EBRS). We use the TSI data set from 
NGDC web site available at http://www.ngdc.noaa.gov/stp and combined 
observational data from National Geophysical Data Center Solar and Terres- 
trial Physics (2012). The importance of UV/EUV influence to TSI variability 
(Active Sun/Quiet Sun) was pointed by Krivova and Solanki (2008). There 
were indicated that up to 63,3 % of TSI variability is produced at wave- 
lengths below 400 nm. Towards activity maxima the number of sunspots 
grows dramatically. But on average the Sun brightens about 0,1%. This is 
due to the increase amounts of bright features, faculae and network elements 
on the solar surface. The total area of the solar surface covered by such fea- 
tures rises more strongly as the cycle progresses than area of dark sunspots. 
Some physics-based models have been developed with using the combined 
proxies describing sunspot darkening (sunspot number or areas) and facular 
brightening (facular areas. Call or Mgll indices), see Frontenla et al. (2004), 
Krivova et al. (2003). 

We also analyzed two activity indices which describe rapid processes 
on the Sun - Flare Index (FI) and monthly Counts of grouped solar flares 
(Monthly flare count). According to Solar-Geophysical Data Reports (2009) 
the term 'grouped' means observations of the same event by different sites 
were lumped together and counted as one. Kleczek (1952) defined the value 
FI = it to quantify the daily flare activity over 24 hours per day. He as- 
sumed that relationship roughly gave the total energy emitted by the flare 
and named it flare index (FI). In this relation i represents the intensity scale 
of importance of the flare and and t the duration of the flare in minutes. In 
this issue we used the monthly averaged FI values. 

The Wolf numbers is a very popular, widely used solar activity index: 
the series of Wolf numbers observations continue more than two hundred 
years. At Figure 2 we illustrated with help of wavelet - analysis the fact that 
the long time series of observations give us the very useful information for 
study of the problem of solar flux cyclicity on long time scales. The result of 
wavelet - analysis (Morley wavelet) of series of observations of average annual 
W is presented in form of many of isolines. For the isoline of the value of the 
wavelet-coefficients are of the same. The maximum values of isolines specify 



the maximum values of wavelet-coefficients, which corresponds to the most 
likely value of the period of the cycle. There are three well-defined cycles 
of activity: - the main cycle of activity is approximately equal to a 10 - 11 
years, 40-50- year cyclicity and 100 to 120-year (ancient) cyclicity. 

3 Activity indices in the cycles 21, 22 and in 
unusual cycle 23 

The recent solar cycle 23 was the outstanding cycle for authentic observed 
data from 1849 year. It lasted 12,7 years and was the longest one for two hun- 
dred years of direct solar observations. This cycle is the second component in 
the 22-year Hale magnetic activity cycle but the 23'''^ cycle was the ffist case 
of modern direct observations (from 1849 to 2008 years) when Gnevyshev- 
Ol's rule was violated: activity indices in cycle 23 had their maximum values 
less then the values in cycle 22 (but according to Gnevyshev-01's rule the 
cycle 23 must dominate), see Figure 2. 

Ishkov (2009) pointed that in this unusual cycle 23 the monthly averages 
values for Wolf numbers during 8 months exceeded 113 and most of sunspot 
groups were less in size, their magnetic fields were less composite and charac- 
terized by the greater lifetime near 2"^^ maximum in comparison with values 
near the 1** maximum.The Wolf numbers reaches its ffist maximum 3.9 yrs 
after the beginning. After the ffist maximum, the index decrease by 14% (of 
that maximum). The two maxima of this index have the same amplitude 

In the cycle 23 the -Fio,7 radio flux and the TSI have the lowest values 
from 2007 to 2009 (the beginning of the cycle 24) all over of these indices 
observation period. The Fioj radio flux index shows the second maximum 
is 8.4% stronger than the ffist one. The TSI (from Earth Radiation Budget 
Satellite) maxima are fainter than those of the other indices because the solar 
irradiance variation of the 23rd solar cycle is just 0.14%. This value seems 
very small but is normal. 

The Flare index - counts of flares has a higher first maximum. This shows 
that the flares can be generated during the ffist maximum, and it seems that 
the generation is decreasing towards the end of the cycle. 

Figures 1 demonstrates that for all activity indices in the 23^^ solar activ- 
ity cycle we can see two maximums separated one from another on 1,5 year 
approximately. We see the similar double-peak structure in cycle 22 but for 
the cycle 21 the double-peak structure is not evident. We see that there are 
displacements in both maximum occurrence time of all these indices in the 
23^^^ solar cycle. 
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Figure 3: Correlation between monthly averages for solar indices versus -Fio,7 
radio flux in the cycles 21 - 23. (a) ^530,3 coronal line flux, (b) Mgll core-to- 
wing ratio, (c) TSI and (d) Counts of flares. 



We assume that the probable reason of such double-peak structures is 
a modulation of the 11-year fluxes variations by both of the quasi-biennial 
and 5,5 year cyclicity. The different time of 1*** and 2""^ maximums appear- 
ance may be caused by the difference in fluxes formation conditions (for our 
indices) at different atmosphere's altitudes of the Sun. 

Figures 1 also shows that for all solar indices in the cycle 23 the relative 



depth of the cavity between two maximums is about 10 — 15%. 

When studied solar activity indices in 21, 22 and 23 solar activity cycles 
we separated out rise phases, cycle's maximum phase, cycle's minimum phase 
and decline phase. In case of linear regression we've found that the maximum 
values of correlation coefficients Kcorr reached for the rise and decline phases 
of the cycles. According to our calculations the highest values of correlation 
coefficients Kcorr we see in connection between W and -Fioj- Correlation 
coefficients Kcorr for linear regression for TSI versus -Fio,? are the minimal of 
all correlation coefficients determined here. Figure 3 also illustrates the high 
level of interconnection between 4 solar activity indices versus -Fio,?- Also we 
see that coefficients of linear regression (slope and intercept) differ among 
themselves for the activity cycles 21 - 23. 

The cyclic variation of fluxes in different spectral ranges and lines at the 
11-year time scale are widely spread phenomenon for F, G and K stars (not 
only for the Sun), see Bruevich & Kononovich (2011). The chromospheric 
activity indices (radiative fluxes at the centers of the H and K emission 
lines of Ca II - 396, 8 and 393, 4 nm respectively) for solar-type stars were 
studied during HK project by Baliunas (1995) at Mount Wilson observational 
program during 45 years, from 1965 to the present time. Authors of the HK 
project supposed that all the solar-type stars with well determined cyclic 
activity about 25% of the time remain in the Maunder minimum conditions. 
Some scientists proposed that the solar activity in future cycle 24 will be 
very low and similar to activity during the Maunder minimum period. Unlike 
this mention Chumak & Matveychuk (2010) predicted the main parameters 
of the new 24*^* cycle closed to the usual activity cycle's parameters not 
similar to the cycle's characteristics during the Maunder minimum. At Figure 
2 for yearly W we also see (if we continue the imaginary line, enveloped 
the maximums) the influence of the century cyclicity to the 24*^ cycle's W 
maximum values. 

4 Changed relation between activity indices 
and Fio,7 in the cycles 21-23 

We analyzed the interconnection between activity indices W, -Fsso.s, Mgll 
core to wing ratio. Flare Index, Counts of flares and TSI versus -Fioj for the 
21**, 22"^ and 23*^ solar cycles. At Figure 4 we presented the interconnection 
between solar indices and radio flux -Fioj in the cycles 21 - 23. 

The linear model corresponds to the linear regression equation: 
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were F^d is the activity index flux, 

Aind is the intercept of a hnear regression, 

Bind is the slope of a linear regression. 



Table 1: Solar activity indices versus Fioj. 
A, B and their standard errors. 



Coefficients of linear regressions: 



Activity indices 


■^ind 


Bind 


Error a a 


Error as 


versus -Fio,7 










W 


-62,28 


1,03 


1,57 


0,011 


-^530,3 


-2,93 


0,084 


0,27 


0,002 


Mgll 


0,258 


0,00013 


0,00026 


0,00000018 


Flare Index 


-7,37 


0,106 


0,51 


0,0036 


Counts of flares 


-269,03 


4,31 


20,14 


0,146 


TSI 


1365,07 


0,0066 


0,044 


0,0003 



In Table 1 we present the coefficients of the linear regressions and their 
standard errors a of intercept and slope values. 

Ishkov (2009) pointed that there was very high level of flared activity 
in cycle 21 and very low level of flared activity in cycle 23. We also con- 
&m this fact (see Figure 3d) that the difference of values Acountsof flares and 
Bcountsof flares determined in 23'''^ cycle from similar values, determined for 
cycles 21 and 22, is more significant among all other indices. The flared 
activity in the 23'''^ cycle almost twice weaker in comparison to 21^** cycle. 

Table 2: Solar activity indices versus -Fio,?- Coefficients of 
polynomial regressions: A, Bl, B2 and their standard errors. 



Activity indices 


A 


Bl 


B2 


error aA 


error asi 


error aB2 


versus Fio,7 














W 


-87,26 


1,45 


-0,0015 


4,86 


0,078 


0,0003 


-^530,3 


-7,38 


0,158 


-0,00026 


0,84 


0,013 


0,000048 


Mgll 


0,25 


0,0002 


-0,0000003 


0,0007 


0,00001 


0,00000004 


Flare Index 


-4,36 


0,057 


0,00017 


1,58 


0,024 


0,00007 


Counts of flares 


-361,4 


5,83 


-0,005 


64,01 


1,01 


0,0035 


TSI 


1364,4 


0,017 


-0,000037 


0,13 


0,002 


0,0000007 



The polynomial model corresponds to the following equation of a second 
order polynomial: 
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Find — ^ind + Blind ' -^10,7 + B2ind ' -^10,7 • (2) 

were Find is the activity index flux, 

Aind is the intercept of a polynomial regression, 

Blind and B2ind are the coefficients of a polynomial regression. 

In Table 2 we present the coefficients of polynomial regressions: A, Bl, 
B2 and their standard errors a of intercept and slope values. 

We have to point out that close interconnection between radiation fluxes 
characterized the energy release from different atmosphere's layers is the 
widespread phenomenon among the stars of late-type spectral classes. Brue- 
vich & Alekseev (2007) confirmed that there exists the close interconnection 
between photospheric and coronal fluxes variations for solar-type stars of 
F, G, K and M spectral classes with widely varying activity of their atmo- 
spheres. It was shown that the summary areas of spots and values of X-ray 
fluxes increase gradually from the sun and HK project stars with the low 
spotted discs to the highly spotted K and M-stars for which Alekseev & Ger- 
shberg (1996) constructed the zonal model of the spots distributed at the 
star's disks . 

We've calculated values Kcorr of linear regression for solar activity indices 
versus Fioj for cycles 21,22 and 23. The values Kcorrit) were determined 
for each moment of time t from Kcorr calculation during the time interval 
t-l,5yr < AT < t + 1, 5yr. 

Figure 5 demonstrates the results of our correlation calculations of these 
solar activity indices versus -Fio,7 {K^orr variations during the cycles 21 - 23). 
We can see that all the K^orr values have the maximum amounts at rise and 
at decline phases. It's easy to estimate the value of period of Kcorr cyclic 
variations as 5,5 years approximately. We assumed that this new cyclicity 
(characterized with period's value equal to half length of 11-year cycle) is 
important for the successful forecasts of the solar activity indices fluxes. 

The cyclic behavior of Kcorr can be explained by following assumption: 
we imagine that some activity index flux depends on time t by the expression: 

FUt) = i^tf ^— '^(t) + AFi^it) . (3) 

were F^^^ aroun ^^^^ ^^ ^j_^^ background flux rising continuously with increas- 
ing of solar activity and AF^^(t) is the additional flux to the overall flux 
from the active regions. 

The previous correlation study allows us to consider that F^^^ s'-o™ ^^^^ 
and AF^^it) are the linear functions of the background and activity regions 
levels of solar activity. In our case we choose the radio flux Fioj as the best 
basal indicator of solar activity levels: 

10 
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poac^grouna^^^ = Oi + 61 ■ Fi''o7^™""''(t) . (4) 

A^nf (t) = 02 + &2 ■ AFf,^{t) . (5) 

The coefficients Oi and bi vary from a2 and 62 in different power for our 
different activity indices. For Wolf numbers this difference is small, but for 
0,1 - 0,8 nm background and count of flares the difference between oi, 61 and 
a2, &2 is more significant than for W and TSI. 

During the rise and decline cycle's phases the dependence Findif) versus 
-^10,7 (^) is approximately linear because coefficients a and h from Table 2 
(which described all the cycle) are close to ai and hi (see Figure 3 and 
Figure 4) and relative addition flux from active regions AF^^f (t) is neglect 
with respect to F^^^ s'-o™ ^^^^ g^ additional flux from active regions cannot 
destroy a balance in the close linear correlation between -Fj„d(t) and -Fio,7(t) 
and respective values of Kcorr reach their maximum from all the cycle. 

During the minimum of activity cycle both values F^^^ g'^oun ^^-^ ^^^ AF^^{t) 
are small, but additional flux from active regions is not neglect in relation 
to background flux that has the minimum values from all the activity cycle. 
Therefore values a and b from Table 2 cannot describe the linear regression 
to a considerable degree in cycle's minimum and values of Kcorr reach their 
minimum values in the cycle. 

During the maximum of activity cycle AF^^{t) often exceeds F^^^ ^'"°"" (t) 
so disbalance in linear regression between activity indices increases and values 
of Kcorr reach their minimum values in the activity cycle too. 

5 Conclusions 

For a long time the scientists were interested in the simulation of pro- 
cesses in the earth's ionosphere and upper atmosphere. It's known that the 
solar radiance at 30,4 nm is very significant for determination of the Earth 
high thermosphere levels heating. Lukyanova & Mursula (2011) showed that 
the for solar 30,4 nm radiance fluxes forecasts (very important for Earth 
thermosphere's heating predictions) there were more prefer to use Mgll 280 
nm observed data unlike usual -Fioj and Wolf number observations. 

The close interconnection between activity indices make possible new 
capabilities in the solar activity indices forecasts. For the successful forecasts 
of maximum values and other parameters of future activity cycles it has been 
required to take into account the century component. At Figure 2 we see the 
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smoothed curved described the maximunis of cycles which shows the century 
cychcity dependence on Wolf numbers cycle's maximums. 

In this paper we found out the cyclic behavior of a calculated during three- 
year interval values of correlation coefficients Kcorr{t) of linear regression for 
W, total solar irradiance, Mgll 280 nm and Counts of flares versus -Fio,? 
during solar activity cycles 21,22 and 23 (see Figure 5). We showed that 
Kcorr{t) have the maximum values at the rise and decline phases - the linear 
connection between indices is more strong in these cases. It means that the 
forecasts of solar indices, based on Fio j observations in our case, will be more 
successful during the rise and decline cycle's phases. 

We also determined that a calculated during three-year interval values 
of correlation coefficients Kcorr{t) are characterized by cyclic variations with 
the period that is equal to half length of period at 11 - year time scale (5,5 
years about). Our study of linear regression between solar indices and Fiqj 
confirms the fact that at minimum and at maximum cycle's phases the non- 
linear state of interconnection between solar activity indices (characterized 
the energy release from different layers of solar atmosphere) increases. 
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Figure 4: Correlation between monthly averages of solar indices versus -Fio,7 
in the cycles 21 - 23. (a) Wolf numbers, (b) ^539,3, (c) Mgll core-to-wing 
ratio, (d) Flare Index, (e) Counts of flares and (f) TSI. 
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Figure 5: Correlation coefficients of linear regression Kcorrif) for (a) TSI, (b) 
Flare Index, (c) Counts of flares and (d) Mgll UV-index versus W and Fiqj. 
Kcorrit) was Calculated during three-year interval. Solar cycles 21-23. 
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